NUMERICAL INTEGRATION FOR HIGH ORDER PYRAMIDAL 

FINITE ELEMENTS 



NILIMA NIGAM AND JOEL PHILLIPS 

Abstract. We examine the elfect of numerical integration on the convergence of high 
order pyramidal finite element methods. Rational functions are indispensable to the con- 
struction of pyramidal interpolants so the conventional treatment of numerical integra- 
tion, which requires that the finite element approximation space is piecewise polynomial, 
cannot be applied. We develop an analysis that allows the finite element approximation 
space to include rational functions and show that despite this complication, conventional 
rules of thumb can still be used to select appropriate quadrature methods on pyramids. 
Along the way, we present a new family of high order pyramidal finite elements for each 
of the spaces of the de Rham complex. 



1. Introduction 

Pyramidal finite elements are used in applications as "glue" in heterogeneous meshes con- 
taining hexahedra, tetrahedra and prisms. Various constructions of high order pyramidal 
elements have been proposed jH [191 El 13 CS] • A useful summary of the approaches 
taken for i/^-conforming elements is given by Bergot et al. [2], who also provide some 
motivating numerical results for the performance of methods based on meshes containing 
pyramidal elements. 

Our aim here is to study the effect of numerical integration on arbitrarily high order 
pyramidal finite elements that approximate each of the spaces of the de Rham complex. If 
they are to be used to implement stable mixed methods, such elements should also satisfy a 
commuting diagram property. One such set of elements was constructed by Zaglmayr based 
on the theory of local exact sequences, [18], and is summarised in [8]. Another construction 
was given by the authors in fT5], and forms the starting point for this work. 

A prototypical (linear) problem is 

(1) For a:V xV and / G V' , find ueV such that: a{u, v) = f{v) ^v^V, 
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where y is a space of functions on a domain, Q C M". One way of obtaining a numerical 
approximation to u is to replace V with some finite dimensional approximation, Vh con- 
structed using finite elements on a mesh whose size is controlled by h. A typical result is 
that the approximate solution converges to the true solution at some rate, 0{h^) where the 
order of convergence, k, depends on the degree of largest complete space of polynomials 
used in the finite element approximation space. 

In general, the bilinear form, a(-, •) and the right hand side /(•) are evaluated using numer- 
ical integration rules. These are additional sources of errors in the approximate solution. 
The theory that describes these errors is now classical and can be found, for example, in 
[Sill]; its objectives, nicely summed up in [5], are 

"to give sufficient conditions on the quadrature scheme which insure that the 
effect of the numerical integration does decrease [the] order of convergence." 

In this paper we will show that the quadratures described as conical product formulae by 
Stroud |16| satisfy the above property for our pyramidal elements. The main challenge 
arises from the fact that the classical theory is only applicable to finite elements with 
approximation spaces consisting purely of polynomials, but pyramidal elements necessarily 
include functions other than polynomials, specifically rational functions (see [15] or [IT]). In 
contrast to the claim in [2], we show that the importance of these functions in constructing 
interpolants means that it is not possible to achieve global estimates of the consistency 
error by summing element-wise estimates that only deal with polynomials. 

Section [2] introduces a framework that will allow us to unify our analysis for discrete 
approximations to each of the spaces of the de Rham complex. We also recall the definitions 
of the approximation spaces for the elements in |15j and the quadrature rules given in [^ . 
In section [3] we show that the conical product formulae are exact for products of all pairs 
of functions from the approximation spaces, including the non-polynomials. The intuition 
from the classical theory would be that this is all that is required. However, in section 
[4] we show that the reasoning behind this intuition is insufficient when functions other 
than polynomials are present. To overcome this, we derive a generalisation of the standard 
Bramble-Hilbert argument. In section [5] we present new families of approximation spaces 
that allow us to take advantage of this generalisation (and which can be used to construct 
new pyramidal finite elements in their own right). Finally, we pull everything together in 
section [6] and show that Stroud's quadrature rules satisfy the desired property for both the 
new and original families of elements. 



2. Definitions 

2.1. Differential forms. In common with Arnold et al. [T] we find that it is natural 
to use tools from differential geometry when discussing approximations to the de Rham 
complex. 
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Let C and define A^*) (O) as the space of differential s- forms on A point, x E has 
coordinates {x^)i=(}..,n. and a given u G A''*)(il) can be expressed in terms of its components, 
u = Uadx°'^ A ■ ■ ■ A dx"" where each Ua G C°°($7) and the multi-indices, a = ai ■ • • 
run over the set, T^, of all increasing sequences, {l...s} — t- {l...n}. 

Define Q^'\n) to be the space of all (covariant) tensors, A : A^'^Q) x A^'^n) C°^{n) 
that can be defined in terms of the pointwise representation, 

(2) A{u, v){x) := A'"^{x)ua{x)vf){x) V«, v G AW(0), 

where we are using the Einstein summation convention, A'^^UaVp := X]a/3eTs ^"'^^q^/?- 
We will insist that is anti-symmetric in the first s and second s components, which 
makes the representation unique. 

A tensor, A G induces a bilinear form on A(^)(r2): 

Jn 

where dx = dx°^^ A • • • A dx"" . 

Let T be a partition of Q where every G T is the image of a simple reference domain, 
K C M", under a diffeomorphism (px ■ K K. On each K, the reference coordinates, 
x = ix^)i^Q „ of any point x E K, are given by x = (j)]^^{x). Given u G A^^\K), the 
reference coordinate system induces a new set of components Ua- Differential forms are 
contravariant, so the components transform as: 

The components of a covariant tensor, A G 9'^*)(fi) transform as: 

Note that (n(x), t;(x))^(^) = A""/^ {x)ua{x)vi3{x) = A"''^ {x)ua{x)v ^{x) is just a 0-form and 
we have the change of variables formula on each element, K: 

(5) {u,v)a,k= / A'^f^UaV0d-x= / A^^UaVB\D(f)K\dk, 

Jk Jk ^ 

where D<pK is the Jacobian of c^k and dx = dx"^^ A ■ ■ ■ A dx"" . 

When n = 2 and n = 3, it is conventional to think of differential forms in terms of proxy 
fields. The spaces A^''''(r2) and A^"'\Q) are always isomorphic to the scalar field, C°°{Q). 
When n = 3, the spaces A(^)(0) and A(^)(0) are isomorphic to the vector field, {C°°{Q))^. 
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For u € A('^)(i7), we denote the components of the proxy field as for i G = {l, . . . , (^)|. 
The isomorphisms for the vector fields are given by 

-""13 • 

\^^3/ \ Ul2 / 

With these identifications, the exterior derivatives, d : A(*)(r2) — t- A(*+-'^)(ri) for s = 0, 1, 2 
become the familiar grad, curl and div. 

As with the differential forms, we will use symbols on the subscripts (and superscripts) of 
proxies to indicate the coordinate system that is being used to determine the components 
of the proxy fields. Given some u G A('')(ri), Uj/ is the ith component of its proxy in the 

coordinate system I. We will also write u' = {uii)i^Xs to indicate all the 

components of the vector (or scalar) field. 

For a coordinate change, x = 4>{x'), the weights appearing in the contravariant and covari- 
ant transformation rules, ([s]) and Q, can be written as the entries of a (^) x (^) matrix, 

. We choose to let to be the weight in the covariant transformation so that, for 

u e A^-'\n) 

(6) ^ (^0 ^) . ., = Vz G Xs. 

The weights can be calculated in terms of the Jacobian, D(/>: 

(7) wf^ = l, wf=Dr^\ wf = \D<I,-^\D4>, njf = \Dcl>-^\. 

The exterior derivative is an intrinsic property of any manifold. This means that it is 
independent of coordinates; equivalently, the exterior derivative commutes with coordinate 
transformation. 

The use of a reference coordinate system is a familiar concept. Shape functions for finite 
elements on simplices are often defined in terms of barycentric coordinates. Using the refer- 
ence coordinate system to examine a shape function thought of as (a proxy to) a differential 
form is equivalent to mapping it to a reference element using pullback mapping. 

2.2. Sobolev spaces. Let the Sobolev semi-norms M^yfc,p(n) and Hj|/fc(Q) = Mvy^.z^j^) have 
their standard meanings. Define semi-norms and norms for any u G A('')(r2) as 

k 

The Sobolev spaces, i?^AW(0) and^(")'''(0) are then defined as the completion of A^'') (fi) 
in the norms ||ti||r,n and ||ii||-^(s),r(f^) = ||M||r-,r2 + M^*||r,r2 respectively. 
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As a short-hand, we will write 'H^'\n) = ^(^)'O(J]). The spaces of proxy fields correspond- 
ing to n^'^^n), n^'^\n) andn^^\n) are the familiar H^{n), if (curl, J]), H{div,n) 
and L2(Q). 

Note0thatir'"+iA(^)(f7) C and in particular = F^+1aW(J1) ^ 

and = H''A(^){n) ^ H''{n). 

For A £ Q^^\0,), we similarly define 

k 

\^\k,oo,Q '■— ^ l^*"* ) ||^||fc,oo,n := l^lr.oo.r2 

ijeXs '-=0 
and define W'''°°e'^''\n) to be the completion of Q^'^Ti) in |M|r,oo,n. 
For a given K and li € A('^)(i^) and A G ©(^^(Er), define the reference semi-norms^ 

Suppose that {Th)h>o ^ family of shape-regular partitions of il, where every K £ Th 
affine equivalent to K and each (px satisfies 

(8) \\D(f)K\\<h and \\D<I)]^^\\ < ^ 

for some p > 1. For any u € 7i^^^'^{K) and A G Ty*'''°°©*^'^)(ir), we have the inequali- 
ties 

(10) ^kh'^"' \^koo,K < \A\,,^,K < Cp'^h'^-'^ 

for some constant C = C{k,n) which is independent of h. These can be deduced from 
the standard scaling argument for Sobolev semi- norms of functions (see, for example, [5j) 
combined with the transformation rules, ([3| and ^ and the observation that ([s]) implies 

that < /i and fg- < f for all 



2.3. Pyramidal elements. From now on we will assume Q C M^. To contain the prolifer- 
ation of indices, we will use the notation (^, r], (^) for the reference coordinates ( 
The reference domain is defined as the pyramid: 

K = {{^,7],O\0<C<l,0<^,rj<C}. 
^When r = 0, this is the observation that H^(nf C i^(curl, Q) and C H{div, Q). 

o 

These are the norms induced by the metric in which the reference coordinates are orthonormal. They 
are used in the scahng argument in section pi 
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We have chosen to restrict our analysis to affine maps (pK so each element of the mesh, 
K E Th, will be a parallelogram-based pyramid. We will refer to a general such K, as an 
affine pyramid. 

As in ipLSj we will also use the infinite pyramid, 

Koo = {{x, y, I < X, y < 1, < z < oo}, 

as a tool to help analyse and define the pyramidal elements. The finite and infinite pyramids 
may be identified using the projective mapping, 



(11) ^-.Koo^K 

(12) (f): {x,y,z)^ 



X y z 
1 + z' 1 + z' iTz J ' 



which can be thought of as a change of coordinates and so, for any element, K, induces 
the infinite pyramid coordinate systerr^ defined as x = (f>~^x. We shall usually write 
X = {x,y,z). 

The corresponding weights in the change of coordinates transformation rule can be calcu- 
lated explicitly: 

(13a) = 1, 

'1 

(13b) w^}^ = D4>-^' = (1 + z) I 1 




(13c) w^^^ = \D(t)-^\D(p= {1 + z) 

(13d) = \D<p-^\ = {l + z)\ 

The approximation spaces for the finite elements presented in |15j are defined on the 
infinite pyramid using k-weighted tensor product polynomials, Q''if^'"'[x,y, z], which are 

tensor product spaces of polynomials, (^''"^'"[x, y, z], multiplied by a weight -r. That 

_ (l + z)'= 

is, Q^'jf^'"' is spanned by the se 

T, < a <l,0 <b <m,0 < c< n} . 



{l + zY 



For each family of elements on the infinite pyramid, an underlying approximation space is 
defined for each order, k > 1. 



so-called because 2; — >■ oo as rj, Q — >■ (0, 0, 1) at the top of the pyramid. 
%f I, m OT n is negative then Q*'™-'" — {0} 
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• i^^-conforming element underlying space: 

(14a) Uf = qI'"''-^ e span ^ 



(1 + ^) 

H(curl)-conforming element underlying space: 

77(1) _ ^fc-l,fc,fc-l ^fc,fe-l,fc-l ^k,k,k~2 
- Vfc+l X Vfc+i X Vfc+l 



• H(div)-conforming element space: 

77(2) _ ^fc,A:-l,A:-: 
L^k — ^k+2 

(14c) ^k-i 



jj{2) _ ^k,k-l,k~2 ^k-l,k,k-2 ^k-l,k-l,k-l 
^k — '^k+2 ^ '^k+2 ^ 




(l + 2)fc+2 

where s{x,y) £ Q''^^'''[x,y], t{x,y) £ Q'''^^^[x,y]. 
• L'^-conforming element underlying space: 
(14d) Z^f =Q^;i''=-i'^-i. 

For an element defined on a pyramid, K, the underlying approximation space, Uj^ (K) is 
defined as the space containing all the s-forms whose components induced by the infinite 
pyramid coordinate system lie in ^^^'^''j^ 

(15) Ut\K) = {nG A(^\K) : (n.)^^^^ G U^^} . 

— — (s + l) 

By inspection, it can be seen that the exterior derivative d -.Uj. Uf^ is well defined, 
and so, since d is independent of coordinates, the exterior derivative on the spaces on each 
element, 

(16) d:Ui'\K)^Ut'\K) 
is also well defined. 

A full explanation of these spaces is provided in [15j. Some motivation may be seen from 
the following lemma. 

^We are using coordinate transformations here, but in [TS], the underlying spaces are defined as the 
puUbacks, ni"\K) = [{(t>~^)*v : v G wi'^j and Ui"\K) = {<j>*iiV : v G lli°\k)}. 
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Lemma 1. For a given K and s € {0, 1, 2, 3} let u & (K). Each component u~- (where 
i £ Is) of u in the reference coordinate system satisfies 

(17) n.o</>GQf'^ 
This means that 

(18) u\^\k) C'H'^'\k). 



Proof. The relationship between the representations of u in the reference and infinite pyra- 
mid coordinate systems is given by equation ([6|): u o (f) = u, where the weights, , 
are given by (13a)-(13d). To estabhsh (17), each s G {0,1,2,3} needs to be dealt with as 
a separate case. 



When s = 0, the weight, w^^^ 
s = 1, inspection of (14b) reveals that U]^"' C Q^+i 



7(0) 



1 and it is clear from (14a) that U]^' C Q^^^^'^ . When 



7(1) 



X Qk+i X Qk+1 ■ The weight, 



'1 

w^}'^ = (1 + z) I 1 



( 1 ) — , .k .k /'~\k.k.k /^k.k.k mi 

so w\ u G Q, X Q, X Q,' . ihe cases s 



2 and 



yX y 1 + Zj 
■s = 3 follow similarly. 

bmce 

Qr C ^°^(^oo) each u- o is bounded on Koa, which means that is bounded 
on K and therefore Ui is bounded on K. Hence ||tt||o,ii' is finite. By (16), du G U^^~^^\k), 
so ||du||o,_ft- is finite too and u G T-L^^\K). □ 



In order to construct pyramidal elements that are compatible with neighbouring tetrahe- 
dral (and hence polynomial) elements, subspaces of the underlying approximation spaces, 

lA}. {K) are identified that contain only those functions whose traces on the triangular faces 
of the pyramid are contained in the trace space of the corresponding tetrahedral elemenfl^l 
These approximation spaces are denoted W^'{K). 

Each approximation space is equipped with degrees of freedom that induce a linear inter- 
polation operator, 

(19) n^^;^ : -?^(^)'i/2+^(i^) ^ ui'\K), e > 0. 

which completes the definition of the finite elements. The necessity of the extra 1/2 + e 
regularity can be seen as a consequence of taking point evaluations at the vertices of the 
pyramid for the s = elements. It is necessary for both the projection based interpolants 
of [7j and the more explicit construction given in [15]. 



The trace spaces for the tetrahedral Lagrange, Nedelec edge and Nedelec face elements are given in 
|13| . By construction, the traces of the underlying spaces on the quadrilateral face of the pyramid already 
match those of the hexahedral elements 
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Given a mesh, 7h, for Q, we can assemble a global approximation space for n^'^n), 

(20) V['^ ={vG n^'^n) : v\k G G Th}. 

The element-wise interpolation operators respect traces on the boundary of the pyramid, 
i.e. 

is) 

tTu\sK = ^ trll^j^nla^ = 0, 

so we can define a bounded global interpolation operator IL^j^^ : '^^'^^'"'^/^'^'^(rj) — t- V^^'' by 
(njf)n)|i^ := ni%{u\K) for ah K G %■ 

2.4. Conical product rule. Quadrature rules on the pyramid can be deduced as special 
cases of the conical product rule presented by Stroud [13 IS] . Stroud defines the quadrature 
scheme for any continuous function, / G C{K), 

(21) S{f) :=^/(ei(l-0),Ci(l-0),^OAiA,W. 

He shows that given n > 0, a sufficient condition for S{p) = Jj^ p dx for any polynomial, 
p G P"'(x), is that the two one-dimensional quadrature schemes given by the points and 
(^i with respective weights Aj and /i; satisfy 

(22) y,K9{C^) = C g{x)dx yg G P", 

(23) V/Ui/i(Ci)= / {i-zfh{z)dz yheP"". 

i -^0 

The A; + 1 point Gauss-Legendre quadrature rule can be used to generate and A, that 



make (22) exact for polynomials of degree 2k + 1. The k + 1 point Gauss- Jacobi scheme 



based on the Jacobi polynomiajj -P^+i^ generates d and /Xj that make (23) exact for 
polynomials of degree 2k + 1. We denote the quadrature scheme for K based on (21 ) that 
uses these points and weights as S^, ^. The error, 

Eu,kU) ■■= S.^^if) - i f{x)dx. 

Jk 

will be zero when / G 



'''The Jacobi polynomials, Pi°''''\s), n >0, are typically defined on the interval [—1, 1]. Under the change 
of variables, s = 2i — 1, they orthogonal with respect to the weight (1 — t)°t^ on the interval [0, 1]. 
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When f E C (K) , where K is a pyramid equipped with a change of coordinates (/)k ■ K — > 
K, we can define the quadrature and error functionals: 




where / = / o (px, i.e. the expression of / in the reference coordinate system, x. 



3. Pyramidal approximation spaces and quadrature 

If u and V are polynomials of degree k then their product, uv G P^'^ so from Stroud's work, 
we know that 5^ x{uv) = /j^ uv. In this section, we shall prove the stronger result: 

Theorem 2. Let K he an affine pyramid; fix k > 1; let s G {0, 1, 2, 3} and let A G Q^{K) 
be a constant tensor field. Then for for any u,v & u'^\k), the quadrature scheme Sk,K 
exactly evaluates the product, (u,v)a,k, i-S- 

Sk,K{A''^UiVj) = {u,v)a,k- 

To do this, we first need to understand exactly which functions our quadrature scheme 
integrates exactly. 

Lemma 3. Suppose that f is a function defined on a pyramid, K , and that the represen- 
tation of f in the infinite pyramid coordinate system: f = f o o lies in the space 

Q2k+i'^''~^^'^''~^^ ■ Then the quadrature scheme, Sk,K is exact for f : 



Proof. It suffices to consider functions p with a representation in the infinite pyramid 
coordinate system: 




P{x,y,z) 




0<a,b,c<2k + l, 
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since these monomials span the space (52^+1'^'^^^'^'^"^^ • In finite reference coordinates, p 
has the form p{^, 77, () = ^°'rf'{l — ^)c-a-6^ ^^^^ using (21 ): 



Sk,K{p) = Sk{\D4)K\p) 

i j I 

= \D(j)K\ [ s^ds [ t^dt / (1 - CT^'^dC 
Jo Jo Jo 

The last step is justified because each of the sums is a quadrature rule applied to a polyno- 



mial of degree <2k + l and so we can apply (22) and (23). Apply the change of variables 
^ = (1 — C)s and 7] = {1 — ()t to obtain: 

Jo Jo Jo 

= \D4>k\ / pi(,,ri,C)dx 
Jk 

= I pdx. □ 
Jk 

We can now prove Theorem [2] where, in fact, we will only need Lemma [3] to be true for 

? ^ ^2k,2k,2k 1 ■ u ■ u r ^2k+l,2k+l,2k+l 

J S (^2^ ' which IS a subspace 01 V2yt+i 

Proof of Theorem^ Let u,v £ ujf\K). A £ Q^'^K) is a constant and so, by the first 
part of Lemma [l| in infinite reference coordinates, the function {u,v)a satisfies: 

Hence by Lemma [Sj 

Sk,K {A'^UjVi) = / A'^UjVi = {u,v)a,k- 
Jk 

□ 



Observe that for the spaces U^^\k), the integrand, A^^UjVj G ^2^-2'^*^ ^'^^ ^' could 
in fact use the scheme Sk-i- 
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4. Numerical integration and convergence 

Let a : Ti^'^n) x Ti^'^n) ^ M be an elliptic bilinear form and let V C n'^'\n) be chosen 
so that the problem of finding u G V such that 

(26) a{u,v) = f{v) WveV 

has a unique solution for any linear functional, f £ V . A discrete version of this problem 
is to find Ufi G Vh such that 

(27) ahiuh,v) = fiv) WveVh, 

where Vh is an approximating subspace of V and ah approximates a using numerical in- 
tegratiorj^ . When Vh is assembled using polynomial elements of degree k, the analysis 
of the effect of the numerical integration is classical; good expositions may be found in 

011!. 

For an example, take an elliptic bilinear form a : Hq{Q,) x Hq{Q,) — )■ M, defined as 

(28) a{u,v) = / A{du,dv) 

Jn 

where A G W^'°°Q^^\Q) and is uniformly positive definite. Assume that Vh C H^{Q) 
is some approximation space assembled using kth order polynomial finite elements and 
that there exists a numerical integration rule, Sh,k,Qi-) which satisfies Sh,k,n{diudjv) = 
J^{diudjv) for any i and j and all pairs of functions u,v G Vh- Let ah{u, v) = Sh,k,n{^idu, dv)). 



It is shown in [Sj page 179] that the solution of (27) will satisfy the error estimate: 

\\u - UhWi < Ch''{\u\i^_^_-^ + \\A\\k,^\\u\\k+i)- 
This result is contingent on an estimate of the consistency error: 

(29) sup r r <Ch \\A\\k,oo\\u\\k+l, 

Wh&Vh \\Wh\\l 

where Hh '■ Hq{Q) — t- Vh is an interpolation operator. The constant C = C{Q,,k) is 
independent of h. 

More generally, an analysis for mixed problems can be found in [9j. The conclusion is the 
same: in order to preserve an 0{h^) approximation error, each bilinear form must satisfy 
an 0{h^) consistency error estimate. 



The key ingredient in the proof of the consistency error estimate, (29) is a local esti- 
mate: 



We choose not to consider the effect of approximating /(■) by some /h(-) using numerical integration 
because it is no different on the pyramid than for other elements. Error estimates may be obtained by 
applying the standard argument and using Theorem [2] 
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Theorem 4 (See [5j, Theorem 4.1.2). Given a simplex, K € Th, assume that for any 
polynomial if) G P'^^~'^{K), the quadrature error, Exitp) = 0. Then there exists a constant 
C independent of K and h such that 



(^W^^'^iK), yp, qGP^{K) 

\0,K 



\EK{A{dp,dq))\ < Ch^\\A\\k,oo,K\\dp\\k-i,K\dq\r 



□ 



This theorem is proved by combining a scahng argument with the following famous result 
from |3]. 

Theorem 5 (Bramble-Hilbert lemma). Let Q C M" be open with Lip schitz- continuous 
boundary. For some integer k > andp G [0, oo] let the linear functional, f : W^~^^'''^{n) — )• 
M have the property that y-ip E P'^{Q), f{'ip) = 0. Then there exists a constant C(0) such 
that 

where |MIvf'=+1'P(c)' operator norm. □ 

In our more general framework, an analogous statement to Theorem [4] might be: 

Proposition 6. Let K e Th be a pyramid. Let s G {0,1,2,3} and A G W''''=^e^''\K). 
Then 

(30) yv,w€Ul'\K) 

(31) \EK,k{A{v,w))\ < Ch^\\A\\k,oo,K\\v\\k-l,K\\w\\o,K 

The problem is that, unlike the situation for purely polynomial spaces, we cannot differen- 
tiate basis functions arbitrarily. We do not have the inclusion, IA^\k) C T-L^^^'^^^lK) for 

k > 3. As an example, take the IA^\k) shape function associated with the base vertex, 
(1,1,0): 

(32) ^(^,r?,C) = ^. 
The L^ norm of its third partial ("-derivative. 



IM) 



(33) Al^ / d^dr,dC 



JO JO 

1 9 



(1-C) 



is infinite. 
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This means that a direct apphcation of the argument in [5l section 4.1] would fail when we 
attempt to use the Bramble-Hilbert lemma (Theorem [5]) to obtain the estimate 

(35) li^^^u ,,<C|n|,^ VrG{0,...,A:}. 



r,K 



An attempt is made to avoid this problem in [2] by using the additional projector vr^ : 
H'+^iK) P' satisfying 

G P''{K) TTrP = p 

on each element, K. This allows element-wise estimates to be established. Unfortunately, 
there is no conforming interpolant onto element-wise polynomials for pyramidal elements 
(see [15] or [I7|). In particular, there will be discontinuities at the element boundaries, 
which means that ||n — vrj.n||i^f7 cannot be bounded. The alternative interpretation of vr^ 
as a global projection onto polynomials would not allow the element- wise estimates to be 
obtained. 

Our solution starts with the observation that not all of the members of each 1J^\k) behave 



as badly as the function v defined in (32). Some are polynomial and others are rational 



but can be differentiated more times before blowing up. For example, we will see in the 



proof of Lemma 13 that u(^, rj, C,)^"^ £ H^^'^{K). So, we start by developing an analogue of 



(35) that, formally speaking, allows us to retain as much regularity as possible. 



Theorem 7. Let C M" he an open set with Lipschitz boundary. Fix a > and let k > a 
be an integer. Suppose that: 

• C -ff"(J7) is a finite dimensional space which includes all polynomials of degree 
k; 

• n : ff"(17) — )• is a bounded linear projection; 

• There exist Vr C H'^{0,) for each r G {0, . . . ,k} such that we can decompose 

i?^ = vb e • • • e Vfc. 

Meaning that for a given u G H^{Vl), the interpolant, liu G R^ , may be decomposed into 
unique functions, Vr G Vr, 

Ilu = vo-\ \-Vk. 

Then we have the following estimates for some of the functions, Vy-: 

• For each r satisfying a <r < k: 

(36) \vrl<C\ul. 

• //, additionally, <zVr, where the space P^ consists of polynomials of homoge- 
neous degree, r, then for each r satisfying a < r + 1 < k: 

(37) \Vr\r < C + \u\^ . 
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Proof. For a given r > a, write Wr = VrU P^^^. Wr C R'^ so we can let '■ ^ Wr be 
any surjective linear projection, is a linear map between finite spaces, so the operator 
(I — ^'roll) : H'^(Q) — )• Wr C H'^{Q) is bounded. Also, since both and IT are projections, 
and P'-^ C P'' C R'' we see that P*" ^ C ker(/ — o H). The Bramble- Hilbert lemma 
gives 

||(/-^'^on)?i||r <c\ui. 

By the definition of Wr, we have (^r ° n)ii = Vr +p for some p G P^'~^ . So |n — -y^ — < 
C|ti|^, which implies 

\Vr\r < \u\^ + + IpI^ = (C + 1) . 

The proof of (37) follows a similar argument. The operator (/ — ^I'j. o H) : H'^~^^{Q) — )• 



Wr C H^{VL) is bounded because r + 1 < /c. The additional condition, P^ C means 
that P'' C Wr and so P'' C ker(/ - ^-r o H). □ 

5. A NEW FAMILY OF PYRAMIDAL APPROXIMATION SPACES 

As identified in |2], the space lA^^ is sub-optimal in that there exist smaller spaces which 
contain the same complete space of polynomials and which are compatible with neigh- 
bouring tetrahedral and hexahedral elements. Here we will identify subspaces, TZ\^\k), 

of each of the original approximation spaces, IA^\k) that can be used to construct finite 
elements with the same approximation and compatibility properties and that still satisfy 
a commuting diagram property. This would be an interesting exercise in its own right but 
within the context of this paper we shall see that the importance of these spaces is that 
they support a decomposition in the manner of Theorem [7] whose components still have 



enough "room" for us to apply a Bramble-Hilbert type argument in Lemma 15 



We start the construction of these spaces in the infinite pyramid coordinate system using 
spaces of A:- weighted polynomials, Q^^'™"', which we define in terms of basis functions (J_^^y 
where a, 6 and c are non-negative integers. 

(38) qL''™^ = span < , ^ , : c<k, a<c + l — k, b<c + m — k 

[{1 + zf 

These spaces can be characterised via a decomposition into spaces of exactly r-weighted 
polynomials, 

k 

[l,m] _ /7^^ ^r+l-k,r+m-kfl 



(39) Qr = ©Q: 



r=0 

It is also helpful to observe that (J_^^^Y ^ £,""q^{l — (^Y~°--^ under the coordinate transfor- 



mation, (r/, ^, Q) = 0(x, y, z) given by ( 11 ). So if the representation in the infinite pyramid 



coordinate system of some polynomial f{x) is / G Q^'™' then / is at most degree k in 
(^, r/, C) and at most degree I and m in (^, C) and (r?, Q) respectively. 
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Now define the spaces TZ^^'^ as 



X gf;*; '1 X {0}) e {Vu : G Qjf'^'^}, 



(40a) 






(40b) 






(40c) 




= ({0} 


(40d) 


nf 


- oi'-' 



Jfc-l,fe-l] 

'k+2 



{vxn:nE(QiV/'''><Ql+V'^><W)} 



The decomposition in the definitions means that exterior derivatives, d : 7^^*^ — )• TZ^^'^^^ are 
well defined. The gradient is injective on Q^^'^^ /M; the curl is injective on {q^^j^I'^^ x Q^^^i ""^^ 

and the divergence is a bijection from ^{0} x {0} x Q^^^2'^ to Q^^j^^'^ ^\ so the se- 
quence, 

R — > nf ^ Tef ^ — ^ 

is exact. The following three lemmas relate these spaces to the Ul ' spaces. To avoid the 
proofs distracting from our main argument we have postponed them to Appendix |A] 

Lemma 8. The spaces TZ^^^ are subspaces of the ti^k^ : 

ni'^'^ui'^ V5G{0,1,2,3}. 



□ 



In fact, for A; > 2 the TZ^i^"^ are strict subsets of the U^k\ 



For a given s G {0, 1,2,3} and k > 0, we define the approximation space on a pyramid, 

(s) 

K, as those differential forms whose infinite coordinate representation lie in 7^^ : 

7^«(K) = {nGA(^)(i^):(n,)G7^^^)}. 

fs) 

These spaces still contain all the polynomials that were shown to be present in the Uj", 
spaces. Specifically: 

Lemma 9. If K is an affine (i.e. parallelogram-based) pyramid then, for k >1, 



□ 



Just as with the original spaces, ujf\K), the new spaces are compatible with Nedelec's 
elements, which were first outlined in |l4j. 

^c.f. the original approximation spaces, ( 15 I 
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Lemma 10. Let K be a pyramid. For each s G {0, 1, 2} there is a trace operator that takes 
elements of 'H^'^\K) to some distribution on the boundary, dK . The image of TZ^^\k) 
under this operator consists of all traces of elements of'H^^\K) whose restriction to each 
triangular or quadrilateral face of K is the trace of a corresponding kth order Lagrange, 
edge and face approximation function on a neighbouring tetrahedron or hexahedron. □ 

Note that the original approximation spaces, ujf\K) were defined by expUcitly identifying 

the subsets of underlying spaces, Uf^ (K) which have polynomial trace spaces. For the 

TZ^j^\k), the polynomial trace property is inherent and this additional step is not required. 
Lemma [TO] is why we need rational functions in our spaces. For example, there is no 
polynomial whose trace is the lowest order bubble on one triangular face and zero on all 
other faces. 

From Lemmas Isl and [Tol we see that: 

Corollary 11. The new approximation spaces are subspaces of the original approximation 
spaces. 

□ 



We can reuse the interpolation operators from the old spaces, (19), to create interpolation 



operators for the new spaces. Since the trace spaces of TZ^/^^ are the same as we 
just need to define projections S^*]^ : U^\k) — TZ^^\k) that do not change the trace 
data. 

For u G IA^\k), define s[[^^n to be the minimum of the functional v i— t- \\d{v — ■u)||o over 
the admissable setj^ 

|w G n^k''^ ■■ u\aK = v\3K and {u - v,dw) = Vu; G n['~^\K)^ . 

Note that is projection-based interpolation of IA^\k) onto 1Z^^\k) and the well- 

posedness of this minimisation problem is established in [7]. 

Now define the maps : H^'^^'^/'^+^K) TZi'\K) as 

V^-^^ ^k,K — —k,K ° ^^k,K- 

(s) is') (s\ is') 

Since both and ^ commute with d, so does j^- In fact, if 11^ ^ is a projection 

is) 

based interpolant, then so is j^- 

^^In the case s = 0, should be interpreted as 
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As with (20), for a given k, we can assemble a global approximation space, 
(42) Sj^^ ={v£ n'{n) : v\K G n^^\K) yK G %} 

and define a global bounded interpolation operator <1>[^'' : ■^'•^^'^/^"'"'^(ri) — )■ S^j^^ by (<l>jj*'*ii)|i^ = 
^^^^k{u\k) for all K G T^. 

We can construct a decomposition for these spaces that we can use with Theorem [7j 

Definition 12. Given a pyramid, K and s G {0, 1, 2, 3} define, for each r > 0, the subspace 
of all the s-forms in TZ^j^\k) whose components are exactly r -weighted when composed with 



(t>:K^^K. 

Note that although the domain of o is Koo , the condition is on the components in the 
reference coordinate system, v--, rather than the infinite pyramid coordinate system vj. In 

effect, what we are saying is that each ?i!j:'^^{K) is spanned by s-forms whose components 
are linear combinations of functions 

(43) e(C,77,C) = eV(l-C)'-'^-'' 

where a, 6 < r + 1. 

Lemma 13. For an affine pyramid, K and for each s G |0, 1, 2, 3} and k>l, each of the 
spaces X^^l {K) satisfy the criterion for Vr from Theor 

Additionally, the semi-norm \ -\^ ^ ^-^ actually a norm on each space X^^^{K). 

Proof. Let u G X^^^{K). Each can be written in terms of functions, e(^, rj, Q) = ^"/^^(l — 
(^y-a-b ^ When a + h > r, these will be rational functions with a singularity aX Q = 1. We 
need to understand their differentiability on the finite pyramid. Let 7 = (71,72,73) be a 
multi-index. The partial derivative. 



7 , In fact. 



e 



dx'y 

where C = C{j, a, b, r) is a (possibly zero) constant dependent only on 7, a, b and r. Hence 

(44) j (^^y = c 1^ 1^''^ l^'U^^^^'''v^'''^H^-Cf^^^^^ 

(45) r(l_^)2(r+l-7i-72-73)^^ 

Jo 

This integral is finite ifr + 1 — I7I > —1/2, soeG H^"^^^^"^^^^ (-^)- By affine equivalence of 
K and k,u£ H^''+^/^-'i (K) C K'+^iK). 
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Finally, (43) shows that each e{^,rj,C,) is either a rational function, or a polynomial of 
degree exactly r, so |e|^^ 7^ 0. Hence \u\^ 7^ and is a semi-norm on X^!^^{K). □ 

Lemma 14. For an affine pyramid, K and for each s G {0, 1, 2, 3} and k > 1, each of the 
spaces TI^^\k) may be decomposed: 



Proof. The decomposition (39) makes the claim look plausible. The details are left to 



Appendix |Xj □ 

6. The effect of numerical integration on the pyramid 

We are ready to assemble all this machinery to prove a version of Theorem |4] for pyramidal 
finite elements. The first step is to establish an error estimate for each of the spaces in the 



decompositions in terms of the reference norms. Recall that in (24) we defined Sk^xi') ^ 
the A;th order quadrature scheme for the pyramid, K, and that we call the error functional 
for this scheme Ek^x{')- We will also use the pointwise representation, A{u,v) = A^^UiVj 
given in ([2]). 

Lemma 15. For any s G {0,1,2,3} and an affine pyramid, K, let v G xj:'^^{K), w G 

n['\K) and A G VF^+i.°°G(^)(^)- Then the error in the evaluation of the bilinear form, 
{v,w)^j^ using the scheme Sk^xi') can be bounded in terms of the reference (semi-)norms 

(46) \EkMA(v,w))\ < C\Dcl)K\\A\^_^^^^^^\vl^^\\w\\^^j^ 

where C = C{k) is a constant that depends only on k. 



Proof. We can transform the error functional onto the reference pyramid using (25). 

(47) Ek^K{A{v, w)) = Ek^K {A'^ViWj) = E^j^ (\D<t^\ A^hiw^^ = \D<I)\ E^ j^ (^A'h.w-^ 

We are able to take \D<j)K\ outside the integral because (px is affine. Define the linear 
functional G G W''-''+^'^e^''> {K)' as 

(48) G{B) = E^j^ (^B'K-w.^ yB G W^~^+^'^Q^'\k). 
Since Sk{-) takes point values of its argument, 

\G{B)\ < C^B^hrW.^^^^ < Cmk-r+l,oo,KMoo,KMoo,K- 

Furthermore, all norms are equivalent on the finite dimensional spaces, X^^^{K) and 



13 



^ ^ is a norm for xj:'^^ . So G is continuous 



TZ^j^\K), and, by the last part of Lemma 
and ||G|| < C'|{'|^^||t()||Qj^. All of the equivalences of norms are done on the reference 
pyramid, so the constant, C depends only on k (in particular, it does not depend on K). 
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From the definition of xj:^^, we know that each v-^ o (p £ Ql-'^^'^^^'^ and by Lemma jlj and 
Corollary 11 w~-o(j)£ Q^^'^ for each j € X^. Now suppose that B is polynomial of degree 



k — r, i.e. each component, B^^ £ P ^ for each i,j £ Ig. Then B^^ o (f) £ Q 
can assemble these facts to see that 



[k—r,k—r\ 



k—r 



. We 



B'^v^wA o </> = iB'^ o {v. o ct>) iw^. o </> £ QllXl 



->2fc+l,2fc+l,2A:+l 



So, by Lemma 



the quadrature error, E 



k,K\B''v.w. 



0. Therefore, P ''' C kerG and 



we can apply Theorem [s] (the Bramble-Hilbert Lemma) to obtain 



\G{A)\<C\A\ 



For some constant C = C{k). Substituting (48) and (47) gives the desired result. 



□ 



We can now apply a scaling argument to get an element-wise estimate on the quadrature 
error. Recall that we defined the interpolation operator, : ?^(^)'i/2+'(K) 'R-k\K) 
in 

Lemma 16. Let K he an affine pyramid satisfying the shape-regularity condition, Q, for 
some p > 1. Fix s £ {0, 1, 2, 3} and take k>2. Then 

(49) Vn £ H''A^'\k),w £ n['\K) and A £ w''+^''^e'-'\K) 

(50) \Ek,K{A{^i%U,w))\ < (C/I'^+I + 0(/l^'+2)) \\A\\k+l,oo,K\W\\k,K\Mo,K 

where C = C{k) a constant dependent only on k. 



Proof. Use the decomposition given in Lemma 14 to write 

^k^K^ = VQ + --- + Vk where Vr G X^ti^). 



By Lemma 15, we know that for each r G {0 . . . fc}, 

(51) \Ek,K{A{Vr, W))\<C \D(I)k\ l^lfe-r+l,oo,J^ \'"r\r,k II^Ho,^' 

The interpolation operator is bounded on ^^'•^^'^/^"'"'^(i^) which is a subset of //3/2+.^{s)(^) 
so Theorem [T] is applicable with a > 3/2. Pick some a £ (3/2,2] so that when r > 2 we 
can use the first estimate, (36), to obtain: 

(52) \Ek,K{A{vr, w))\<C \D(I)k\ l^lfc_,+i,oo,i^ l^lr.i^ ll^llo,K- 

Now apply the inequalities ^ and (10) to the semi-norms (and norm) on the right-hand 
side to obtain 

hr+s us 

\Eu,K{A{vr,w))\<C\D^K\h^-'-^^-^'p''\A\^ 



\k—r+l,oo,K I ^ , \l/2 



r,K ,^ , ,1/2 



K\ 



Ch^^^ |^|fc_r.+l,oo,i^ l^lr.i^ \\wh,K, 

where the generic constant, C still depends only on k. 
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When r = 1, we can similarly apply the second estimate from Theorem [T] given in (37) to 
obtain: 



(53) \Ek,K{Mvi,w))\<Ch''+'\A\,^^^^j,(\u\,^j, + h\u\^^j,) \\w\\o,K. 
For r = 0, note that H^^ollci- < ^ll'"ll3/2+£,x - ^ {\'^\o,k + \'^\i,K + I'"l2,ii')' 

(54) \Ek,K{Aivo, w))\ < Ch^+^ |^|fc+i,oo,i^ {\<,K + h \u\,^K + \u\2,k) \Mo,k 



so 



Summing over the Vr, we obtain (50). 



□ 



Summing these errors over each element gives an estimate for the global consistency error 
due to the numerical integration (we shall ignore the 0{h^~^'^) terms). Recall that in (42) 

we defined the global approximation space, sl"^ cn^'\n). 

Theorem 17. Let s € {0, 1, 2, 3}, k >2 and assume that is constructed using a shape 
regular mesh, % and finite elements, TI^^\k) for each K £ %■ Let A £ Ty''+^'°°G('*)(r2) 



an 



d u £ H^^^'^{Q). Then the interpolant ^\^^u G S'^^ satisifies 



sup 



< /i'^^"'"||^||A;+l,oo,Q||'u||A:,Q 



\Wh\\0 



Where we define {v,w)A,h,k,n ■= Hk^Th ^^^^ (^(^'"^))- 



Proof. Let Wh e s\^\ 



KeTh 



< Ch'^^^ ^ ||^IU+i,oo,_ft:||^i||fc,_fs'||'w^/i||o,K 

KdTh 

1/2 



< Ch^+HA 



fc+l,oo,f! 




Dividing through by Htf/iHon gives the result. 



^ \\wh\\l,K 
□ 



In the proof of Lemma 15, the important condition for w was that w~-o(f)£ Q^^'^ . So, by 



Lemma 



we could equally well have taken w G ujf\K). Furthermore, S'f^'' C V^*'' means 



(s) (s) 

that <I>^ n G . Hence we have a consistency error estimate for the global approximation 

(s) 

spaces based on the original elements: 
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Corollary 18. Under the same assumptions as Theorem 11 , let be constructed using 



finite elements based on the approximation spaces, U^\k). Then the interpolant ^)^'u 
satisifies 



sup 



k+i,oo,n\\u\\k,n- 



□ 



The error estimate may be applied to more general bilinear forms because of the commu- 

(s) (s+1) 

tativity d o 11^ = 11^ o d. For example, the consistency error for the elliptic bilinear 
form, (26), is 

{d^''l^\,dv)A,n - {d^f^u,dv)AXk,Q. 



sup 



a{^fu,v)-ah{^';;'u,v. 
Iblli 



(0) 



< sup 



(0) 



< sup 



\dv\ 



{^''j^^du,w)A,n - {^h'du,w)A,h,k,n 



(1) 



\w\\o 



< Ch^^^\\A\\k+i,oo,n\\au\\k,n 

< Ch^^^\\A\\k+i^Qo,n\\u\\k+i,n- 

A final note: as with the classical theory, the error estimates decay like 0{h^^^) but these 
are emphatically not /ip-estimates. The degree, k enters into the constants in several places, 
which is to be expected from arguments that rely on the Bramble-Hilbert Lemma. 



7. Conclusion 

The conventional finite element wisdom is that a A;th order method requires a A;th order 
quadrature scheme. We have shown that this is still true for some high order pyramidal 
finite elements, but that the non-polynomial nature of pyramidal elements requires some 
unconventional reasoning to justify the wisdom. 

In the process, we have demonstrated new descriptions of families of high order finite 
elements for the de Rham complex that satisfy an exact sequence property. We will examine 
these elements in more detail in future work, but a couple of notes are worth recording 
here. 

• The approximation spaces for the first family in the sequence, TZ^^^ (K) are the same 
as Zaglmayr's elements, as described in [8], and which [2j describes as optimal with 
respect to their dimension and compatibility with neighbouring elements. 
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Lemma 



shows that the TI^^\k) spaces contain polynomials corresponding to the 
tetrahedron of the first type. Zaglmayr has constructed pyramidal elements con- 
taining polynomials corresponding to both types of tetrahedron, but only those 
corresponding to the second type are presented in pj. It would, clearly, be inter- 
esting to compare our TZ^^^ {K) spaces with the construction for the first type. 



A. Properties of the new approximation spaces, TZ^ 
In this appendix, we have collected proofs of various Lemmas in Section [5} 



Proof of Lemma [S[ The inclusions 

(Jy^^ Q[/,m] J— / Ql,m,uiin{l,m} — 1 _|_ QO,0,min{i,m} J Qi,m,min{i,m} 



can be verified from the definition, (38). By the first inclusion, Q^f''^' ^ Qfj''^'^ + Q'I'^''', 



which gives the s = case: TZ^i^^ C U^^^ . 

The s = result implies VTZf^ C VUf\ Thus, since VU^^^ C U^k\ we have VQ^'''^ C U'j}\ 
which establishes the result for the second space in the decomposition for Tz2^\ given in 



(40b). To deal with the first space in this decomposition, apply (55) and the definition of 



Ui^ given in (14b), to obtain 



Qfc''^ X QK"' X {0}) ^ (Qt'i'^'-' - Qlil'^'-' X {0}) cUi'\ 



(2) 

The ,5 = 2 case may be established similarly. The space TZi is defined via a decomposition 



"(2) 

ofU)^ by taking curls of the s = 1 result. The first space is dealt with by applying (55) 
directly to the definitions. 

Another application of (55) gives TZ^j^^ = Q^^_^^''^ ^' ^ ^ '^ ^ =^^k ^- 



into two spaces, (40c). The second space in this decomposition can be seen to be a subset 



□ 



Proof of Lemma^ Since P^ is preserved by affine transformation, we can work in the 
reference coordinate system, x. Recall the components of the proxy representation of some 
u G A^'^\K) in this coordinate system are denoted u-^, where i G Ig- We will need to show 



that if all the components, u-- £ P^ (or, for s = 1,2,3, ^) then u G TZ^^\k). This is 
equivalent to showing u G TZ^^\ which we will do using the transformation rule, ([g]), along 
with the explicit weights associated with the coordinate change (j) '■ ^oo K given in 
(TSab-pMl). 
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We start with the case s = 0. Let u G A^^\K) be any polynomial, u(^, r/, () = ^'^rf'{l — CY 
where a + b + c < k. Then 



u 



w 



(0) 



-1 



u o 



^ayb 



1 



a+b+c ^ 



n 



(0) 



Similarly, for s = 3, take u G A(3)(i^) as u{C,rjX) = ^"rf'il - CY ioi a + h + c < k - 1. 
Then 



(1 + ^) 



7^ 



(3) 



The s = 1 case involves a little more work. Let u G A^^^ {K) have polynomial components, 
lij G P^^^ . We can find g G A('^)(i^) with representation q G P*^ such that v = u — Vq has 

third component (in reference coordinates), = 0. By the result for s = 0, g G TZ^^\k), 
and so (by pOb] )) Vg G t6^\k). We need to show that v G T6^\k). Both and t;^ are 
in P^~^ . Suppose first that = ^"•rf'{l — CY where m'.= a + b + c<k — 1 and = 0. 



w 



(1) 



-1 



V O (j) 




a+l„,b-l + 



rV 



m+1 _j_ ' 



(1 + 2)™ + ^ 



Compare this last expression with the definition, (40b]), to determine that v G Note 
that when a = m (which includes the case a = /c — 1), the first term vanishes, because 
6 = and 1 — = o|^ An identical calculation establishes the same result when v | = 

and ^2 = rr?^(l - C)''- 

For s = 2, the change of coordinates formula for u G hY^\K) is 



(56) 



w 



(2) 



U O 



{l + zf 




Suppose that = ^"r/ (1 — ^)'^ with m 

t 



contribution to u is 



0,0 . Let p 



a + 6 + c<A; — 1. Apply (56) to see that the 

1 x'^y'' 

+ 2 (l + 2)™ + 2 



^ ^ ^ G and observe that 



9p 

'92 



and 



— — — — - — G q£_2^''^ (^^^ case b = m implies that a = 



dx 



m+2 (l+2)™+2 



In other words, (5", 0,0)* is an exact 1-form. 
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and therefore ^ = 0, so the final inequahty in (38) is not violated). Hence 



w 



(2) 



1 - c)^ 













1 o (/) = V X 


P 


l-l 










\dx/ 





Polynomials in the second component can be dealt with similarly. When = ^""q {1 — CY, 
the contribution to u is yjp^^y^, (i+^)m+3 , (i+^)m+2 j • Hence 



(2) 



Note that 








V X 



1 

m + 2 









{H-2)'"+2 1 








H 






v 


J 








1 _ a+b+2 \ x'^y^ 
^ m+2 j (l+z)™+2^ 



{1+zy 



G Q|f ""^'^ unless a = m or b = m, but in these cases, (l — 







□ 



(2) 



Proof of Lemma 10, This can be proved in an identical manner to Lemma [ref] in [15] for 



the original spaces, h(jf\K). We will just give a sketch here. First we need to show that 

the restrictions of the traces of the TZ^i^\k) functions to each face lie in the trace spaces 
of the corresponding tetrahedral or hexahedal approximation spaces. Secondly we need to 
show that any valid trace can be achieved by some member of TZ^^\k). 

Convenient definitions of the tetrahedral and hexahedral spaces may be found in [13j . As 
an illustration, observe that members of the TZ^^"^ which are non-zero on the face y = of 
the infinite pyramid can be expressed in terms of monomials p^^, where a + c < k, which 

map to which will span all polynomials of degree k on the face r/ = of the finite 

pyramid, which is precisely the trace space of the fcth order Lagrange tetrahedron.. 

The second step is proved by demonstrating a linearly independent set of shape functions 
with non-zero traces that is the same size as the set of external degrees of freedom. This 
can be done by showing that it is possible to achieve the lowest order bubble on each face, 
edge and vertex of the pyramid. The sets of shape functions presented for the IA^\k) in 
[H] also suffice for the n\!\K). □ 



Proof of Lemma 14' Each X^^^ is a subset of T^)^' , so 

AfgWe---eAfW(K)c7e«(K) 

For the reverse inclusion, we will deal with each s G {0,1,2,3}, in turn. For every s G 
{0, 1, 2, 3}, the transformation rule, Q, gives uo (p = w^^^u. 

^^Just as earth-shattering, this is the observation that (0,0, 5")' and (0,0, 77'')* are exact 2-forms. 
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For 0-forms, the weight in the change of coordinates formula w^'^ is equal to 1 so any 



u E 7^^**^ (K) satisfies u < 



uGTZ 



(0) 



The decomposition, (39) gives 



which is a subset of Qq'^'^ 



■■Q 



fc+l,fc+l,0 



0,0,0 



Q 



k,k,0 



For the cases s = 1 and s = 2, we will consider a basis for TZ^j^\k) and show that each 
element, u, of the basis is a member of xj'^j}{K) for some r G {0 . . . k}, which amounts to 
showing that each w- o (j) ^ Qr^"^'^^^'^ ■ 



From the definition given in (40c) it's natural to consider three cases for an element of a 

From the 



basis for t6^\k). First suppose that u € (^(5^+/'"' x x Oj with -uj = 

\k—l k] 

definition of Q)^j^i ' we see that < a < c — 2 and < 6 < c — 1 and so 2 < c < /c + 1. 



Jfc-i,fc] 



Then w^j^'^u 



x"-y" 



T,0, 



ij and so each u-- G Q'^^l''''^ C Qr where r = c — 1 G 
{1 . . .k}. The second case is when u G (O x Q^^]^i x j and the reasoning is identical to 



r,r,0 



the first. Finally suppose that u = Vp where p = G Qj^ 

Vp = 0. So we can take c > 1 and see that each entry of 



[k,k] 



When c = 0, p = 1 and 



(1)- 

w, u 



a 



r,.a-lyb 



"(l+z)<=-i 



V(Q + fe-c) (i+4c-i / 



is in Ql-^^'^'^^'^ for some r G {0 . . . A;}. 

When u G TZ^^\k), lets start with the case u G ^0 x x Q^^_^_2'^ SL,nd write Ug 
^ ^ Again, it is simple to check that each of the entries in the vector w^^^u 



(i+zy 



is in 



r+l,r+l,0 



for some r G {0...k}. Now suppose 



that u = V X V where v G [Q^^_^_i'''^ x x 0) with 



(l+z)- 



From the s = 1 case, we 



know that c > 2 and so its straightforward to verify that each of the entries in 



(2)~ 



(1+^) 




(1+2)"= 



' + = 

-cx°-y'' , -bx^y'' 
(l+z)':-2 (l+z^-^ 



are in Qj-"''"'^'^'''^''^ for some r G {0 . . . A;}. The argument for u = \/xv with v G x Q^+i 
is the same. 



[fc,fc-i] 



X 
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Finally, u G TZ^j^\k) means that u G Q^^^^''^~^^- The weight w^^^ = (jq^ so u o cj) = 
(i+z)^ ^ ^ Qt-i'''~^^ the reasoning is the same as the 0-form case. □ 
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